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The Coulomb gas partition sum is rewritten in terms of the path integrals 

formalism. It is shown that perturbation theories based on the Mayer expan- 
sion and on the path integrals method lead to the identical results. The well 
known Debye-Hukkel result for the case of 3D Coulomb plasma is completely 
rederived. An analogous result is obtained for the case of the Coulomb gas 
with dipoles. This result can be considered as a generalization of the Debye- 
Hukkel approximation. Other possible generalizations including the range 
of interaction potentials for which the Debye-Hukkel approximation can be 
applyed are discussed. 

I. INTRODUCTION 

The problem of calculating of the Gibbs partition sum is well known since the pioneering 
works in classical satatistical mechanics. The simplest example of precise calculaton of the 
Gibbs free energy is the classical ideal gas. The more general problem of classical statistical 
mechanics is the evaluation of partition sum for a nonideal system. The main difficulty is 
the calculation of the configuration integral 



This problem has no explicit solution. There are several approaches to analyse the ther- 
modynamics of non-ideal systems with short range interactions: Mayer's method, correlation 
function method, integral equation method, renormalization group method e.t.c. 
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Some additional difficulties appear while analysing Coulomb-like systems because of the 
long range type of the interection. Since the pioneering work of Debye, the standard solution 
of this problem is the Debye-Hukkel (DH) mean field approximation |1|,0] where the Debye- 
Hukkel screened interaction appears instead of Coulomb one. This method is widely used 
in different fields of physics to analyze 3D plasma-like systems. 

Howewer, the extention of DH theory to the case of more complex systems is still a 
challenge to the statistical mechanics. On our opinion, the main reason for this is the 
absence of the consistent point of view on DH approximation. Really, a number of ways to 
obtain DH result are known. First of all this is the linearization of the mean field equations 
which is the simplest and well known way to obtain DH screening. Another way to obtain 
the DH result is a cycle approximation |l|] applied to the Mayer's diagrams. Howewer, like 
other chain approximations in statistical mechanics, this method leaves open the question 
how to choose the relevant chain class. We will discuss another method based on the special 
property of the configuration integral and the analogy with multiple Gaussian integrals. 

Despite the fact that the first attempt to apply method discussed below to continuous 
systems was undertaken in 1959 0, to our knowlege this approach has never been used as a 
regular method to analyse classical systems. At the same time some exotic relations between 
different correlation functions were obtained by using this approach (See for example @]). 
The analogous method has been used to prove the equivalence of a 2D Coulomb gas model 
to the Euclidean Sine-Gordon model 0. This equivalence was considered as a tool allowing 
to obtain some new results in quantum field theory from thermodynamics properties of the 
2D Coulomb gas. 

There exists some papers in statistical mechanics [^], where authors try to obtain some 
new results for Coulomb systems by using the Sine-Gordon analogy and the technique of 
the quantum field theory. But due to some difficulties some well known results for the 3D 
plasma were not completely rederived in these papers. 

Some purely mathematical papers published in 80-th should be specially mentioned [0. 
In these papers the 3D Coulomb system has been analysed on the base of its equivalence to 
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the 3D Sine-Gordon model. A strict mathematical theory was developed for this case. As 
it was shown, the Debye-Hukkel (DH) result can be easily reproduced in the lowest order 
of this approach. This is the simplest and the most consistent method to obtain the well 
known DH result. 

In those papers, however, the emphasis was made on the mathematical problems, and 
it is difficult to see the physical sense of the tranformations. Also there are no discussions 
on possible extension of this method to the case of more complex non-Coulombic systems. 
The present paper is an attempt to overcome some of these difficulties and to present the 
application of Gaussian integral method to continuous systems. 

In the section |I| the basis of the Gaussian integral formalism applied to the models 
of classical statistical mechanics is described. Section |ITT| is dedicated to the perturbation 
theory, that is fully analogous to the usual quantum field theory's one. The correspondence 
of this theory to the usual Mayer's expansion in powers of the activity z is analysed in detail. 
Some mathematical restrictionscaused by the features of our method are also considered. 
As it is shown, these restrictions are strong enough. This explains, why the DH-like mean 
fied theory can not be applied in the case of non-Coulomb systems. 

In the section |V| the expansion in powers of density is presented. It is shown that the 
lowest order of this expansion gives precisely the DH approximation. The exact reproduction 
of this well known result permits us to expect that the method used in this paper can lead 
us to some new results in the case of more complexes systems. One such example is given 
in section 0. The results of the present paper are summarized in 0. Some advantages and 
defects of the method are reviewed in the conclusion. 

II. GAUSSIAN FORMALISM FOR THE CONTINUUM CASE 

Let us consider the functional: 

y" ^X(^)exp|-i^ y" xir)drH{r ~ r')dr'x{r') + ip J p{r)xir)dr \ = I[p], (1) 
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J\f = I Dx{T)expl—— I xi^')drH(r — r')dr'x(r') 



where: 

Af = [ D\(r\er^n[ - ^ 

- is a normalizing factor. H{r) is a kernel of some linear operator. Dx{r) is a mesure on 
the space of real functions. In other words, in Eq.(|l]) we integrate over all pathes x(^) as 
one usually do in quantum field theory and quantum mechanics. 
Now we can introduce a function V{r) satisfying the condition: 

V{r - r')dr'H{r' - r") = 6{r - r"). (2) 

This condition means that V is the inverse operator relative to H. (Here and later in this 
paper we suppose H and V to be translationally invariant) 

With the help of V{r — r') we can transform (P to the form: 



exp |— /5/2 J p(r)drV{r — r')dr' p(r')^ = I[p{r)] 



'^exp'j/5/2 / x{'^)drH{r — r')dr'x{r') + 1(3 I p{r)x{r)dr\ . (3) 



N 

Later we shall examine the conditions under which the above transformation is valid. 

The main aim of this section is to associate Eq.(|^) with the statistical mechanics of a 
classical system. For this purpose we consider a system consisting of = + charged 
point-like particles. The microscopic charge density of the system can be written as: 

P,{r) = Y.^{r- r,) + - r,). (4) 

We also suppose that the interparticle interaction is central - SiSjV{rij) (where Sj, Sj - are 
the charges of particles), so that we get for the full energy: 

^ i^j=l ^ i^j=l i=l i=l 

pq{r)V{r - r)pq{r')drdr' - {N+ + N^)V{0) (5) 

Using Eq. (^ and Eq. (0), we can write Gibbs factor exp(— /3f/(ri . . .tat)) of the system 
through the path integrals: 



/ JV_,_ JV_ JV_|_ JV_ 



Here we denoted for simplicity: 

(X, H,x) = J x{r)drH{r - r')(ir'x(r'), 
(X,P<7) = / X{r)pq{r)dr 



Equation will be the main for our later examinations in this paper. Applying analytical 
techniques of path integrals calculations to the Eq. we hope to receive mathematically 
based results for thermodynamics of the system in question. Talking about mathematically 
based results we mean the results received directly from the basic physical formulas, us- 
ing usual mathematical ideas without additional physical supposition. Our purpose is to 
demonstrate the main way to obtain physical results, while the details of a mathematical 
consideration is the topic of interest for mathematicians. 
Let us note the following: 

a) Interparticle potentials V{r) ordinary contain sigularities at r — 0. Moreover, there 
is a large and very important class of potentials weakly decreasing at r ^ oo. (Coulomb 
potential in ci = 2, 3 is the simplest example) Due to these effects divergences can appear 
in H{r) and its integrals. (See Eq. (^)) This may cause each component in the right hand 
side of Eq.(|]) to be infinite. On this stage of consideration we shall ignore this fact by using 
for our purposes the most general notation. For mathematical completeness we can suppose 
all divergences caused by the form of V{r) to be somehow regularized. The regularization 
parameters should be taken equal to at the final stage of calculations after taking the usual 
limit A^-^oo,f2^oo,p = const. We suppose that the final result does not depend on the 
type of regularization. 

b) Let us consider the range of the applicability of Eq. (^) . Despite a seeming generality, 
this formula has some very important restrictions. To understand them we can write the 
discrete matrix analogue of Eq. (||) : 
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- 2 y ^ xiA.ijXj-\-i y ^ ^k^k 

lE^-A^'^j / dxi...dxNe 



e =^ ^ . (7) 



(ixi . . . (ixTve 
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It is wel known that this equation is correct only if is a positively definite form. Oth- 

erwise integrals in the numerator and the denominator of Eq. are diverging. Therefore 
Eq. is correct only for positively definite operators H{r — r'). Due to the translational 
invariance of H, we can write this condition with the help of Fourier transforms as: 



H{q) = J dre''^''H{r) > 0; 



Viq) = ^ > 0. (8) 
H{q) 

We suppose here that the interparticle interaction V{r — r') = V{\r — r'\) is central, so that 

V{q) = V{-q) = V*{q) 

Coulomb potential is the simplest example satisfying (H): 

Viq) = ^, rf = 2,3; (9) 

Yukava (screened Coulomb) potential is the second example: 

V{q) = ^^, d = 2,3; (10) 
+ A"^ 

Another (slightly exotic) example is the intervortex interaction in a thin superconducting 
film i: 

V(q)= , d = 2 (11) 

27rAg(g + l/A)' ^ ^ 

Here (po = hc/2e is the fiux quantum, A is a penetration depth in two dimensions. Statistical 



mechanics of vortices interacting via this potential is analysed in [pi)[[r2 . 

c) One should pay the attention to the factor e^^*-^-' in Eq. (|^). As far as we know, the 
most authors ignore this factor. To exclude it the potential 

nr) = (12) 

^ ^ I 0, r < a ^ ^ 



was used by Edwards instead of Coulomb potential {d = 3). It is easy to show that 
Fourier transform of Eq. (|I2|) is not positively definite, so that the correctness of is 
doubtful. 

The most strict method is proposed in mathematical papers (see for example where 
the interaction: 

1 e"''/^ 

Vnir) = -- e- (13) 

r r 

is used, instead of Coulomb potential. The final result should be received after taking the 
limit £ — > 0. This provides V{q) to be positive, and V{r) equal to zero as r — >• 0, so that 
our factor ~ e^'''-"-* is equal to 1. 

It is easy to see that Eq. ( pTSD is no more than one of different ways to regularize V{r). 
We can consider a more general method to smooth singularities in V{r): 

VR{ri — = J dr'dr"n{ri — r')V{r' — r")n{r" — r2), 

with Vr(0) 7^ 0, so that the factor e^^^^'^'^f'^ cannot be omited. We shall see later that this 
coefficient is very important in the mean field and perturbation theories, based on Eq.(^. 

III. DIAGRAMMATIC EXPANSIONS 

To apply the Gauss-analogy (|^) to the problems of statistical mechanics let us consider a 
system consisting of "positively and negatively charged" particles, having absolute activities 
Zj^ and The term "charged particles" is employed here in unusual sense: The interparticle 
interaction may be not Coulomb, but it depends on some special charge-like characteristics 
of particles. For example the interaction between two vortices in a supercondoctor depends 
on mutual direction of their magnetic fields. 

The grand partition function for such a "plasma" at a temperature 1/ j3 can be written 

as: 




Using Eq.(^ (or Eq. (j^)) we can rewrite the above equation: 



N+ N- 
Z, Z_ 



(-)- 



,l3NViO)/2 



-(3/2ix,H,x)^i(3 J pg(r)x(r)dr f^-^^-^ 



We will suppose our system to be "neutral", so that z+ = = z. Keeping in mind Eq.( 
we get: 



c 



/ (^-P/'^ix,H,x) + 2z J cos px{r)dr 



(16) 



where z = ze'^^^'^^^'^. For the case of Coulomb gas, where V{q) ~ ^, and H(r — r') ~ 
5{r — r')V^ the equation (|1^) becomes similar to the Euclidean Sine-Gordon theory. This 
equivalence is a well known fact for more than two decades . To calculate approximately 
C we can use a perturbation theory in the spirit of quantum field theory |^ . For this purpose 
we need to consider again: 



Z\J] = 
and to define: 



at" 



exp f — j J{r)V{r — r')J{r')drdr' 



(17) 



< f{x{r)) >H= I ^/(x(r))e-^/^(^-^'>^). f{jj^^)Z[J] 



(18) 



j=o 



In this case we obtain for (: 
( = (exp (25 / cos f3xi''")dr 



H 



oo 1 

E 



n=l 



2k 



dr 



Z[J] 



■ (19) 



j=o 



Following methods widely used in quantum field theory, we represent each term of the 
sum in (|19D as a diagram with n vertex connected by k lines. It is easy to show that an even 
number of lines are closed in each vertex. Each vertex (i) contributes an integration over 
dvi, the factor —i3V{rij) corresponds to each line connecting vertices i and j. 

Quantum field theorists know that ( (Eq. (|19|)) is infinite in the limit Q ^ oo. The 
infiniteness of ( (Eq. (|I4|)) is also known from statistical physics. To our good both cases 
need only the value of \n(/Q. This permits us to exclude all disconnected diagrams, like 
one usually acts in quantum field theory. 
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Now it's time to examine diagrams, containing tadpoles. Tadpole is a line {—pV^ru)) 
which begins and ends at the same vertex, so that the resulting muliplying factor is propor- 
tional to V{0). The term "tadpole" is well known by field theorists, but is used very seldom 
in classical statistical mechanics. 

Let us consider a vertex with / tadpoles and j non-tadpoles lines outgoing from it (exter- 
nal lines with respect to vertex). The combinatorial prefactor in the integral, corresponding 
to the diagram with such a tadpole is: 

Ml 

/!j!2'' 

where 2k = 21 + j is a full number of lines outgoing from the vertex, (note that j must be 
even) 

Summing by a number of tadpoles / of the vertex with fixed number j of external lines, 
we can prove, that all tadpoles of the vertex can be absorbed into a factor e~^^^'^^ bound to 
this vertex. This is equivalent to the renormalization of z to its original value: 

z = e^^(o)/2z ^ e-^^(°)/2~ ^ ^_ (20) 

Now the role of the factor e^^^(°) in Eq. (§) becomes clear. The existence of this factor 
enables us to compensate the singular contributions of tadpoles in diagrams. To compensate 
these tadpoles in quantum field theory the normal ordering of the initial Gamiltonian is 
usually used (see for example [|ll|). In the case of statistical mechanics we have no base 
neither to use normal ordering, nor to suppose our initial parameters to be infinite, so that 
the singular factor appearing in Eq. (|^) is the only mathematically correct way to avoid 
singularities bound with tadpoles. 

At the same time it is important to note that such a re-renormalization ( [20| ) takes place 
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only if we consider the full interaction cos /3v^.[| In the case of substitution: 

cos(3^= 1 ^ + ^P' 

the maximum number of lines in the vertex becomes equal to 4, and the maximum number 
of tadpoles is limited, so that the factor e^^^^^ can not be fully compensated by them. 

After the discussion of tadpole-compensating it is not difficult to show by analogy, that 
summing of multiple lines, connecting two vertexes allows us to change each set of diagrams 
with multiple lines by one diagram with "bold" line, where bold lines are assosiated with 
Mayer's function fij = e"^^*^'"'^^ — 1. Here we shall not concentrate on the proof of this fact. 
The analogous result is described in 

The main conclusion from the above consideration is that the perturbation theory derived 
from the path integral representation of the Grand partition sum (see Eq. (|16])) - path 
integral perturbation theory, is equivalent to the Mayer's series derived directly from Eq. 



(p!^), and can be written through the same parameters: z, fij. 

At the same time, despite the similarity of these two evaluation, there are a difference 
between them: really, Mayer's diagrammatic expansion of ([l^) contains two kind of ver- 
tices, corresponding to both (positive and negative) kinds of particles. Factors and z_ 
correspond to these vertices. Each pair of vertices can be connected by two kind of lines, 
depending on their mutiple charges: e~^'^+'^^^ and e~^^** To understand why /+_-type 
lines do not appear in the path integral perturbation theory, we can expand Mayer's function 
fij = e~^'^^^ — 1, in powers of Vij. The difference between /_|_+ and will be seen only 
in diagrams, containing one or more vertices with even number of lines (even vertices). For 
each diagram containing even vertex we can find another graph, differ only in a "sign" of 
this even vertex. Due to the neutrality of our system (2;+ = z_ - see above) these diagrams 
will compensate with each other, giving no contibution to the final free energy. Thus we 



^Here the quantum field theory terminology is used: (x, x) - \s a. Hamiltonian's kinetic term, 
while {cos{(3(p)) is the interaction. 
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have shown that for Mayer's developement of the grand partition sum, after representation 
of fij as series (see for example |jl|]) diagrams with even vertices do not contribute to a final 
result. Summing the rest of diagrams we get the resut identical to the one obtained from 
the path integral perturbation theory. 

In this section we have shown schematically the similarity and moreover, the equivalence 
of both perturbation theories. This result gives us no new methods to calculate thermody- 
namic properties of the system, but is important for understanding different approximations 
obtained from Eq.(|l^) or directly from Eq. (^^). This is not a purpose of this section 
to perform a full consideration of different perturbation theories connected with Mayer's 
series. Here we consider only those results that can by obtained from the path integral 
representation of the partition sum. 



IV. THE EXPANSION IN POWERS OF DENSITY 



To obtain the expansion in powers of density, it is convenient to consider the partition 
sum instead of the grand partition sum, used in the previous section. For the system of two 
sorts of particles with densities p+ and p_ the partition sum can be written as: 



NJ. NJ 



J {dr} exp I -i^ ^ Vij {vij) I 



(21) 



Here {dr} denote the integration over all + particle positions in d- dimensions, A 

,1/2 . 



\~2^) ^ Boltzmann wavelength. We obtain the free energy F: 



-(3F = \iiZ = N^ 



where In Q is the nonideal part of the free energy: 



















. \P+J 









+ lnQ, 



\nQ{N+, N., P,n) = In J 



dvi dr^ 



-I3U 



(22) 



Q is the (i-dimensional volume of the system. 
Now we rewrite exp(— /3f/) using Eq. (0): 
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Q = eM^Vm J £^exp{-P/2{x,Hx)}e''\ (23) 

= iV+ In / ^e*x{r+) + iv_ In / ^g-^^^''-) 

To calculate the path interal Dx we need to diagonalize the quadratic form in the exponent 
of Eqs. For this purpose it is convenient to use Fourier transform as a change of 



variables in path interals |14| 



X{r) = I ^X,e-^^^ X, = I rfrx(r)e^''^ (24) 



x{r) is real function, its Fourier transform satisfies the requirements x* = x~q- 

For our further consideration we suppose x to oscillate around % = 0, so that 5*1 can be 
expanded just to the second order in x- In the case of neutrality (p+ = P-) we get: 

P+ + P-i^ 12 P+ + P- f dq ^ ^ 



^1 = -^\Xo\ - J j^M . (25) 

Path integration over Dxij) should be replaced by: 

DXq = dxo n dReXqdlmXq, 

where g > denotes the production through a "half-space of q. Here it is important to 
remind the mathematical definition of the path integral as a limit of multiple integral over 
function values on the discrets set of points. Keeping in mind the simmetry of l^(r),and 
H{r) we can write: 



J x{r)drH{r - r')dr'x{r') = J j^JXq\'H{q). (26) 
In this case for the exponent in Eg . (^31) we get: 



- /?/2 E ^(^ + ^qWq) - E + 4') - f3/2^Hml- (27) 

Here we use the summation instead of the integration over g-space following the rule: 

dq 1 ^ 
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The notation J2 means a summation over "half -space V of the g-space , so that for each 

q>0 

q in V, —q is not in V. Considering this sum we specially marked a point q = 0. We also 
introduced for convenience: 



Xo = Xo- 



Now we can represent the configuration integral Q using the above formulas: 



Q = exp (^V^(O)) / n d^,d^,dXo 



X 



X exp 



^ E {-PH{q) - (3\p^ + p„)) + ID - [^m)xl 



Analogously, the normalizing factor M can be written as: 
A/" = / n diqdlqdxo X 



X exp 



After simple calculations, we obtain for In Q: 



lnQ = ^r(0) + ^ln. 



mq) 



q>0 

At the final stage the sums should be replaced back by the integrals: 

In Q = - 2 y ln(l + /3(P+ + P-)l^(g)) + j ^ 

The last term in Eq. (|32D reflects the fact that: 

dq 



(27r)' 



lV{q). 



Finally, for the nonideal part of the free energy we obtain: 



I3AF _ 1 

n ~ 2 



^ {ln(l + /3(p+ + p_)V{q)) - Pip^ + p^)V{q)} 



(28) 



(29) 



(30) 



(31) 



(32) 



(33) 



For the case of the neutral 3D plasma (y), this equation is known as a Debye-Hukkel 
approximation for the free energy. The way, presented above ia a new method to ob- 
tain the well known result directly from the Gibbs partition sum without special physical 
suppositions. 
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The fact that Eg. ([33|) is identical to the result of the summation of cycles diagrams 
can be simply explained. Really, the expansion of 5*1 in ( p3D just to the second order of x 
corresponds to the retaining for later summation only those diagrams which have no more 
than two lines in each vertex. The vertices without lines appears in disconnected diagrams 
only, graphs with one line vertexes compensate themselves if the system is neutral. The only 
diagrams leaving for the summation are cycles, and their summation leads to Eq. (^3|). 

At the same time, the path integral method described above explains us, why the result 
(^) cannot be applied to the non-Coulomb like systems. For the Van-der-Vaals interparticle 
interaction V{r), the Fourier transform Vq is usually not positively defined (see section ||). 
Due to this fact we cannot apply all the above considerations to the non-Coulomb-like 
systems. To obtain the analogue of Eq. (|33D in a general case, we need to invent the 
analogue of the Gauss representation (|^) with a sign-alternating quadratic form Aij and to 
extend it to the continuous case. The aim of this paper is not to realize this plan, but to 
show some features of the path integral method in statistical mechanics, and to obtain some 
new results for the simple case of Coulomb gas. 



V. SYSTEM WITH DIPOLAR PAIRS 

The interesting chance appears while using the path integral representation of the Gibbs 
partition sum: we can include clusters of particles in our consideration. For simplicity let us 
consider a system, consisting of positively charged and negatively charged particles. 
A part of these particles is organized in clusters of known form. The intercluster interac- 
tion and their self energy can be calculated by summing "each with each" interparticles 
interactions. This fact can be easily used to include clusters of particles in Eq.(|D. 

As a simple example we consider a system of N^ + N_ positive and negative free particles. 
We will suppose also that there are Nc dipoles. Here a dipole is considered as a pair of 
particles (+, — ) with distance between them to be fixed and equal to a. Such a dipol has 
2d — l degrees of freedom. We can write the microscopic charge density Pq of the system as: 
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N+ N- Nd 

P, = Y: Sir - rt) -Y.S{r- rj) + Y.{5{t - rf +) - 5{t - rf)), (34) 

i=l j=l k=l 

where r+ and r„ are positions of the free particles and r-^* = ± a/2 are positions of 
positive and negative charges of the dipole. 
Gibbs partition sum for the system is: 



N+\ NJ Nd\ J [ n j \ \ ^0 J 



Here z± are activities of free particles and zd - activity of dipoles|. We integrate in ( pSD over 
positions of all free particles (dr^), and over positions of dipol's centers (dr^). In addition 
in Eq. (|35D the integration over all orientations of dipols (^^^ where ooq is a normalizing 
factor is performed. For the case = 2 we have i^ = ^, forrf = 3- ^^ = siR^^. 
Taking into account the "each-with-each" interaction of particles, we write for U: 

U= J2 s,s,V{r,,), (36) 

so that U includes not only ions-dipoles and dipoles-dipoles interactions, but interaction of 
particles inside the dipoles. This fact is important when we define the activity of the dipoles. 
Due to Eq.(|^) we represent through the path integral: 

^3 = exp I \ j exp(-/3/2(x, Hx) + '^2) ^ , ^ • (37) 

Where N = N+ + N_ + 2Nd, and 



^2 = iV+ In / ^e'Mr) + iv_ In / ^e-'^r) + Nn\n [ ±^e'(^^^^'^+-m-^irn~a/2)) .gg) 

J \l J \l J LUq 



Now we expand 5*2 just to the second order of as done in section |IV|. For convinience 
we use the Fourier transform of x{r)- Using: 

x{rD + a/2) - xirD - a/2) = -2i j -^x{q)e~''^''^ sin(qa/2). 



^Here we consider dipoles as the third kind of particles 
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we get for 5*2: 



^2 



(27r)^ 



(39) 



where a function: 



(40) 



is used. In Eq.(|^) the integration is over all directions u of vector a. 

Using the same notations as in section |I^ we can perform the same sequence of calcula- 
tions and for the exponent (similar to Eq. (|27|) ) we can write: 



m E + mil) - P'^^^ E + el) - ^P'PD^ + e|)r(g)- 



g>0 



2 



n 



q>Q 



(41) 



The configuration integral can be written as: 



Q^exp(fr(0)j^ 



^7/11 d^,di,dxo exp U E - P\p+ + P-)- (42) 



Here - is the same normalizing factor as in section ||. Finally we obtain for In Q: 

1 



-/3F, = lnQ 



ln(l + 4/5p^r(0)r(0))- 



(43) 



q1 

2 



dq 
(2^ 



{ln(l + P{p+ + p_)\/, + 4/?pz?r(g)) - (p+ + p_ + 2p^)fiV:,/2} 



where the first term can be neglected as it remains finite in the limit Q ^ 00 

Now we can review additional suppositions introducted in this section while obtaining 
Eq. (|43|). The only special supposition in comparison with section ^ is the definition of 
dipoles and their including in the partition sum (pS]). Other derivations, mathematically 
are fully equivalent to ones made in section Moreover, it is obvious that in the limit 
p/5 Eq.(^) reduces to the nonideal part of the free energy of pure two-components 
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plasma Eq.(^). This fact makes us to hope that the result ( ^31) can be considered as a 
generalization of the Debye-Hukkel result for a two compomemt plasma with dipoles. 

The including of dipole pairs is not the only way to get some new analytical results using 
the path integral formulation in statistical mechanics. For example to describe a hexatic 
phase it may be useful to introduce hexagonal clusters, so that the charge density becomes: 

N+ N- Nh 

p(r) = ^ 5(r - r+) - ^ 5(r - r") + J2 Pah{r - r^), (44) 

1=1 j=l k=l 

Where pah - is a local charge density of the hexagonal cluster of size a: 

6 

Pah{r - Tk) = So6{r - r^) + ^ SnS{r - + a^). 

n=l 

(sj) are charges of particles of clusters considered. 

In this case we shall get a function T{q) that is different from Eg. (^0]) . This function will 
be closely related to the symmetric properties of clusters. 

Analogousely we can introduce a partition sum including many different types of clusters. 
The resulting free energy can be considered as a function of their densities. We may hope 
to obtain the physical properties of our system by minimizing this free energy. For example, 



the free energy Fd (Eq. (^31)) can be considered as a function of free particles (p+ ~ p_) and 
dipoles (pd) densities. If the minimum of such an energy is placed in the region {po ^ p+), 
we can conclude that under these conditions free particles couple into dipole pairs. 

VI. CONCLUSION 

In the present paper a new point of view on the classical Gibbs partition sum is pre- 
sented. This method is based on the representation of the configuration integral through the 
Gaussian path integration. Obtained results precisely reproduce the Debye-Hukkel result 
for the symmetric 3D Coulomb plasma. As it is shown, this method and its extensions 
can be applied to the case of Coulomb and quasi-Coulomb systems. Presented method is 
applyed to the very interesting case of symmetric Coulomb plasma with dipoles and the 
genralization of DH result was obtained. 
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Presented point of view on the statistical mechanics has some advantages and defects. 
The main defect of this point of view is a formahty of computations. Really, all trans- 
formations presented above are simple extensions of the well known techniques of multiple 
Gaussian integrals computation. By the way we were forced to ommit the analysis of some 
questions concerning the divergence of intermediate results. 

Talking about the advantages of the described method, it should be mentioned that we 
got a method allowing us to find the range of the applicability of the DH result (see section 
HI). Also we can see possibles systematic ways to obtain higher order correction to the 
known results. For example it is obvious that the next order of the DH result is the well 
known corrections to the Gaussian approximations of the integral: techniques of calculation 
of these corrections is well developped, and we can try to generalize it to the continuum case. 
Analogously, we can extend DH (or DH-like) result for the case of more complex systems 
with Coulomb interparticle interaction, as is done in sec. 0. 

Method presented in this paper should be extended and analyzed from differents points 
of view: First of all it should be interesting to obtain numerical results using new formulas 
derived above and to compare these calculations with others approximations. We hope 
to get new results by applying the path integral representation to Coulomb systems with 
differents types of clusters,as it is done with dipols in section 0. As it was pointed above, 
the generalization of the Gaussian path integration to the case of non positively defined 
Fourier transforms is needed to widely test the presented method. Another topic of interest 
is the connection of the path integral representation with the renormgroup approach and 
critical phenomena. There is a good chance that the presented non-trivial representation 
of the Gibbs partition sum can give us new approaches and methods to physics of critical 
phenomena. 
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